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We investigated the ferroelectric domain wall propagation in epitaxial Pb(Zr,Ti)03 thin films 
over a wide temperature range (3 - 300 K). We measured the domain wall velocity under various 
electric fields and found that the velocity data is strongly nonlinear with electric fields, especially 
at low temperature. We found that, as one of surface growth problems, our domain wall velocity 
data from ferroelectric epitaxial film could be classified into the creep, depinning, and flow regimes 
due to competition between disorder and elasticity. The measured values of velocity and dynamical 
exponents indicate that the ferroelectric domain walls in the epitaxial films are fractal and pinned 
by a disorder-induced local field. 

PACS numbers: 05.45.-a,47.15.G-,64.60.F-,68.35.Rh,77.80.Dj 
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The physics of surface growth in disordered media 
with quenched defects is of crucial importance to un- 
derstand numerous intriguing natural phenomena [l[ , in- 
cluding contact lines in wetting, surface of epitaxially 
grown films, and magnetic domain walls. In such media, 
elastic forces tend to keep surfaces flat, while defects lo- 
cally promote the wandering, as schematically displayed 
by Fig. 1(a). The competition between elastic and pin- 
ning forces leads to a complicated energy landscape with 
many local minima, which affects the surface growth dy- 
namics under an external force. Recently, there have 
been extensive reports to adapt the fractal concepts to 
surface growth dynamics Q. 

Ferroelectric (FE) domains have been studied for 
past decades because of scientific importance in micro- 
scopic aspects such as multi-domain formation, stability, 
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FIG. 1: (color online) (a) A schematic diagram of domain 
wall propagation in disordered medium. Elastic forces come 
from the curvature of domain wall, and defects work as strong 
pinning sites, (b) Theoretical prediction on the domain wall 
velocity v vs. electric field E in system governed by competi- 
tion between disorder and elasticity effects. Eco represents a 
threshold E. 



and pattern at equilibrium as well as technological appli- 
cability in multi-functional devices such as FE random 
access memories, actuators, and sensors 0, [3]. Quite 
recently, lots of piezoresponse force microscope (PFM) 
studies have provided us microscopic aspects of FE do- 
mains, including inhomogeneous nucleation process Q 
and the fractal nature of their rough surfaces Q. Note 
that most works on FE domains have been focused on 
their static properties. In spite of its scientific and tech- 
nological importance, we have limited understandings on 
how the FE domain wall propagates in terms of time. 

We suggest to prospect FE domain wall from the view 
of nonlinear responses, which follow the predictions of 
the statistical physics on surface growth. Then, the FE 
domain wall velocity v should have a nonlinear behavior, 
shown in Fig. 1(b), under E. At zero temperature T, the 
domain wall remains strongly pinned by local disorders 
until E reaches a threshold value Eco- When E Si Eco, 
it experiences a pinning-depinning transition and starts 
to move with a nonzero velocity v, as represented with 
the red dashed line. Under this depinning regime, 



(E - Eco) 6 



(1) 



with a velocity exponent 9. Under the flow regime, when 
E ^> Eco, v ~ E. On the other hand, for finite T, the 
pinning-depinning transition becomes relatively smooth, 
as represented with the green solid line. Under the low 
E (i.e., <C Eco) creep regime, the domain wall motion 
becomes very slow and can be described by propagation 
between pinning sites due to thermal activation. Then, 



v ~ exp[-{U/k B T){Eco/EY} 



(2) 
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where U is an energy barrier and \x is a dynamical expo- 
nent. Critical exponents of domain dynamics, including 
fi and 9, can identify the universality class and provide in- 
formation on the pinning forces and the fractal nature of 
the rough FE domain walls. Although measurements of 
these critical exponent values for the FE systems are par- 
ticularly important, there are little experimental works 
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FIG. 2: (color online) (a) P-E hysteresis loops at different T. 
(b) Time-dependent normalized switched polarization, Ap(t), 
behaviors (solid symbols) from switching current and (open 
symbols) from PFM measurements. The solid lines show the 
fitting results using the KAI model, (c) Scanned images of 
t-dependent domain growth at E = 0.3 MV/cm. (d) Plot of 
1/to vs. v. This relationship is linear, suggesting that 1/to 
from switching current studies can be used to parameterize v. 



on these values . 

In this Letter, we report our studies on the T- and 
independent nonlinear responses of FE domain wall dy- 
namics in epitaxial Pb(Zr,Ti)0 3 (PZT) thin film. To 
widen the accessible region of T and E, we used switch- 
ing current measurements, combined with direct v data 
determined by PFM images. We found that v follows the 
nonlinear dynamic response, as described in Fig. 1(b). 
We also could obtain values of the two critical exponents, 
[i and 6, from the data in the creep and the depinning 
regimes, respectively. This work provides us new insights 
on how domain walls propagate inside epitaxial FE thin 
films. 

We fabricated 100 nm-thick epitaxial PZT thin film 



on SrRuOa/SrTiOa substrate using pulsed laser deposi- 
tion Q • X-ray diffraction studies confirmed that a high- 
quality, (OOl)-oriented PZT film was grown epitaxially. 
To fabricate PZT capacitors, we patterned the sputtered 
Pt top electrodes with a typical area of 7.5 x 10 3 /im 2 . 
Our epitaxial film revealed a high dielectric stability, suit- 
able to our measurements at high electric fields. 

Figure 2(a) shows T-dependent polarization-electric 
field (P-E) hysteresis curves for a PZT capacitor, mea- 
sured between 3 and 300 K. The saturation and remnant 
P values were nearly constant over a wide T range. The 
systematic T variation of the P-E hysteresis curves comes 
mostly from T-dependent change in coercive field Ec ■ At 
3 K, Eq ~ 1 MV/cm. As T increased, Ec decreased sig- 
nificantly. 

One of the difficulties in performing the dynamic stud- 
ies on FE domains is to measure reliable values of v un- 
der uniform E. Recently, we developed a modified PFM 
technique for a FE thin film with a top metal electrode 
[E 111 • By combining PFM with switching current mea- 
surements, we were able to track wall motions of domains. 
After applying one positive poling pulse (10 V, 50 fxs) to 
pole P, we switched P with a series of negative pulses. 
We then measured the PFM images after all the negative 
pulses. We assumed that the PFM image obtained after 
the negative pulses would be nearly the same as that ob- 
tained after a single pulse, with the width being equal to 
the sum of all the negative pulses Q. Typically, an im- 
age acquisition process using the PFM set-up takes a few 
minutes, which is relatively long in relation to the time 
scale corresponding to the width of the E— pulse used in 
these experiments. To check its validity, we calculated 
the amount of normalized switched polarization Ap from 
PFM images with an area of 5 x 5 /im 2 . These are shown 
as the open symbols in Fig. 2(b). We also determined 
Ap independently from switching current measurements, 
which are shown as the solid symbols. The agreement 
between these data indicates that we could reliably use 
the modified PFM techniques to study domain wall mo- 
tion in our FE film. 

Figure 2(c) shows typical PFM images of time- 
dependent growth for a particular isolated domain, which 
were obtained at room T with E = 0.3 MV/cm. As time 
t increased, the domain size increased. Eventually, it be- 
gan to merge with other domains. From the area of the 
isolated domain, we determined its mean radius. By di- 
viding the change in mean radius with the change in t, 
we obtained the v value. We repeated these procedures 
for more than 20 isolated domains and finally obtained 
an average value for the given T and E. Note that our 
modified PFM study can provide a way to measure v di- 
rectly. However, obtaining sufficient data to plot all v vs. 
E curves would have been too laborious. In addition, it 
is rather difficult to use in the low T region. Therefore, 
it is highly desirable to find another method to obtain v 
reliably in a wide range of T and E. 

The switching current in a FE thin film should be orig- 
inated from the P reversal, whose behavior is governed 
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FIG. 3: (color online) Ap(t) obtained via switching current 
measurements under various E values at (a) 300 K and (b) 
3 K. The solid lines show the fitting results using the KAI 
model. 



by domain wall dynamics. The t-dependent changes in 
Ap for epitaxial FE thin films 9] have been explained us- 
ing the Kolmogorov-Avrami-Ishibashi (KAI) model [Iol |. 
According to the classical statistical theory on nucleation 
and unrestricted domain growth, 



Ap(t) = 1 - exp[-(t/t ) n ], 



(3) 



where n and to are a geometric dimension and a 
characteristic switching time for the domain growth, 
respectively. In the two simplest cases, analytical 
relationships between to and v can be easily obtained. 
When the nuclei of opposite polarity are generated 
at a constant rate under E, t ~ (l/v)^" 1 )/™ [nj. 
Conversely, when all the nuclei are generated instanta- 
neously, to ~ 1/v 1 1 01 ] - Our previous studies on epitaxial 
PZT films showed that nucleation rate is approximately 
proportional to \/t which is much closer to the latter 
case, so 1/to might be nearly proportional to v. 

We investigated the relationship between v and 1/to 
using the switching current response. Our experimental 
Ap(t) data were fitted using the KAI model, i.e. Eq. 
(3), as displayed by the solid lines in Fig. 2(b). This 
prediction was closely correlated with the Ap(t) data. 
We compared 1/to values with v values, which were 
directly measured by PFM. As shown in Fig. 2(d), 
1/to is linearly proportional to v with a small offset, 
indicating that the switching current response could 
provide reliable values of v. 

Figure 3 shows Ap(t) data with varying values of E. 
The solid symbols and lines indicate the experimental 
data and the fitting results using Eq. (3), respectively. 
At 300 K, a sudden change in Ap(t) with E = 0.1 
MV/cm occurred between 0.1 and 1.0 s: namely, at 
to ~ 2 x 10 _1 s. As E increased, a sudden change in 
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FIG. 4: (color online) (a) (Solid symbols) T-dependent curves 
of 1/to vs. E. The dotted and solid lines are guidelines for 
eye. (b) log(to) vs. 1/E curves. The linear fitting indicates 
a dynamic exponent (x ~ 1. The inset shows that UEco/ks 
is nearly independent of T. (c) Plot of log(l/to) vs. log(E — 
Eco) with the Eco ~ 1 MV/cm. From the slope of the plot, 
we found a velocity exponent 6 ~ 0.7. 



Ap(t) occurred over a shorter timescale, i.e., a smaller- 
value of t . Conversely, at 3 K, little changes were 
observed in Ap(t) when E was below 0.8 MV/cm. This 
implied that the FE domain was pinned by defects and 
could not move below a threshold value of E. However, 
Ap(t) starts to change abruptly around E — 1 MV/cm. 
This corresponded to the pinning-depinning transition. 
When E increased, to rapidly decreased. The Ap(t) 
data for 3 and 300 K were similar at E = 1.1 MV/cm; 
t - 5.4 x 10~ 6 and 3.3 x 10~ 6 s at 3 and 300 K, 
respectively. These similar to values demonstrated that 
the T-dependence of v became insignificant under very 
high E region. 

Figure 4(a) shows experimental 1/to vs. E curves. 
Note that these 1/to — E plots resemble those in Fig. 
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1(b). In the low E region, 1/to was strongly dependent 
on T, consistent with predictions for the thermally 
activated creep regime. However, in the high E region, 
the values of 1/to started to merge, indicating a crossover 
to the flow regime. The pinning-depinning transition at 
3 K occurred at approximately 1 MV/cm, which was 
close to the Eq value obtained from P — E hysteresis 
curve in Fig. 2(a). 

From the experimental 1/to data, we obtained 
the value of the dynamical exponent fi for the creep 
regime. From Eq.(2), log(t ) should be proportional 
to {U/k B T){E C o/Ey. The log{t ) vs. 1/E curves are 
shown in Fig. 4(b). The experimental data at a given 
T falls approximately into a linear line, indicating that 
fi is close to 1.0. Considering experimental errors, we 
found fi = 0.9 ± 0.1. This fj, value agrees with previously 
reported values for epitaxial and polycrystalline PZT 
films [111, [Hj]. The inset shows the value of UEco/ksT 
calculated from lines of best fit for several values 
of T. The straight line indicates that the value of 
UEco/ks is nearly independent of T and is about 300 
K-MV/cm, close to previously reported value of around 
400 K-MV/cm from local domain switching data using 
PFM at room T [PH. 

The value of \i reflects the nature of the pinning 
potential in our PZT films. Under a pinning potential 
with a short range (the so-called random bond), onc- 
and two- dimensional domain walls should have the p, 
values of 0.25 and 0.5, respectively [HI]. However, under 
a pinning potential with a long range (th e random field), 
/i = 1.0 regardless of dimensionality pj|]. Recently, two 
conflicting fi values were reported for epitaxial PZT 
thin films; 1.0 and 0.5-0.6 [13, 03. Our study confirms 
that \x = 0.9 ± 0.1 over wide range of T. This \x value 
suggests that the defects in our PZT thin film induce a 
long-ranged local field and pin FE domain walls [lj| . 

From our experimental 1/to data, we also obtained 
the value of the velocity exponent 9 near the pinning- 
depinning transition. By Eq.(l), log(l/to) should be 
proportional to 9 ■ log(E — Eco)- As shown in Fig. 4(c), 
this is approximately the case when T = 3 K. The slope 
of the line for best fit was 9 « 0.71 ± 0.05. We could not 
find any earlier studies with which compare this 9 value 
for FE thin films. 

The #-value should reflect the dimensionality D of the 
surface for an elastic object in a disordered medium, as 
9 = (5 + D)/9 [li|. Using this relationship, we found 



that D w 1.4 ± 0.4 for our PZT thin film. Recently, 
several studies have reported non-integer dimensionality 
of local and static domain walls in FE thin films 0, Il6j |. 
By measuring the local switching of PZT-BiFeC>3 sol-gel 
thin films with liquid PFM, Rodriguez et al. reported 
that D « 1.5±0.1 [6j. By measuring morphology and 
scaling of the domains in epitaxial BiFeOs thin films, 
Catalan et al. reported that D w 1.5 ± 0.1 [16]. Our 
measured values of 9 from the dynamic responses is 
consistent with these D values from the static measure- 
ments. In addition, our findings indicate that the newly 
observed fractal dimensionality of FE domains should 
be originated from local quenched defects. 

We want to point out future studies and possible im- 
plications of our works: (i) in order to fully understand 
the domain wall dynamics in epitaxial FE thin films, 
we need to measure other critical exponents, related 
to the divergence of the correlation length, the local 
variance of the domain wall position, and so on [2j; (ii) 
domain dynamic responses in systems with different 
structural conditions such as polycrystalline FE thin 
films and ultrathin films should be investigated and 
compared with those in epitaxial films; (iii) the complete 
understanding on FE domain dynamics is necessary for 
numerous practical applications, including optimizing 
operation speed of miniaturized FE devices. 

In summary, we investigated domain wall motions of 
epitaxial PZT film over a wide range of temperature 
and applied electric field. We found that the motions 
were likely to be governed by the nonlinear dynamics of 
surface growth in a disordered medium with quenched 
defects. We determined two critical exponents for 
domain wall propagation dynamics, which indicate the 
random field nature of the defects and fractal nature 
of domain walls. Our works provide us a new future 
direction of studies on the domain wall motions in 
ferroelectric materials. 
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